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Background Mutations in KRAS, BRAF, PIK3CA and PTEN expression have been in focus to predict the effect of epidermal 
growth factor receptor-blocking therapy in colorectal cancer (CRC). Here, information on these four aberrations was collected and 
combined to a Quadruple index and used to evaluate the prognostic role of these factors in CRC. 

Patients We analysed the mutation status in KRAS, BRAF and PIK3CA and PTEN expression in two separate CRC cohorts, 
Northern Sweden Health Disease Study (NSHDS; n = 197) and Colorectal Cancer in Umea Study (CRUMS; n = 414). A Quadruple 
index was created, where Quadruple index positivity specifies cases with any aberration in KRAS, BRAF, PIK3CA or PTEN 
expression. 

Results Quadruple index positive tumours had a worse prognosis, significant in the NSHDS but not in the CRUMS cohort 
(NSHDS; P = 0.003 and CRUMS; P = 0.230) in univariate analyses but significance was lost in multivariate analyses. When analysing 
each gene separately, only BRAF was of prognostic significance in the NSHDS cohort (multivariate HR 2.00, 95% CI: 1 .16-3.43) and 
KRAS was of prognostic significance in the CRUMS cohort (multivariate HR 1.48, 95% CI: 1.02-2.16). Aberrations in PIK3CA and 
PTEN did not add significant prognostic information. 

Conclusions Our results suggest that establishment of molecular subgroups based on KRAS and BRAF mutation status is 
important and should be considered in future prognostic studies in CRC. 



Colorectal cancer (CRC) is one of the most common causes of 
cancer- related deaths in the western world (Jemal et al, 2008). 
Distant metastases represent the greatest threat to patient survival 
and about 40% of the patients will die from a metastatic disease. 
Surgical resection is today the basis for curative therapy, but a 
detailed understanding of the biological processes that regulate the 
establishment and progression of a malignant tumour may lead to 
improvements in non- surgical antitumour therapy. Two develop- 
mental pathways of sporadic CRC have been identified: chromo- 
somal instability (or microsatellite stable, MSS) and microsatellite 
instability (MSI). Microsatellite stable tumours are considered to 
arise by copy number gains of oncogenes and loss of tumour 
suppressors, due to numerous chromosomal translocations 
(Grady, 2004). In contrast, MSI tumours show loss of expression 
of mismatch repair genes. They are less often associated 



with lymph node metastasis and distant spread, and MSI patients 
have a better prognosis than stage-matched MSS patients 
(Gryfe et al, 2000; Kohonen-Corish et al, 2005; Popat et al, 
2005; Wright et al, 2005). Additionally, MSI tumours have been 
associated with CpG island methylator phenotype (CIMP) (Ahuja 
et al, 1997), where the groups CIMP-high, CIMP-low or CIMP- 
negative are based on promoter methylation frequency. 
We and others have reported a poorer prognosis for CRC 
patients with CIMP-high or CIMP-low tumours, compared with 
CIMP-negative tumours, especially in combination with MSS 
(Van Rijnsoever et al, 2003; Ward et al, 2003; Samowitz et al, 
2005; Ogino et al, 2007; Shen et al, 2007; Barault et al, 2008; 
Dahlin et al, 2010). 

Signalling through receptor tyrosine kinases in response 
to cytokines, growth factors and hormones is important for 
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maintaining the metabolism, proliferation, survival and motility 
of a cell (Haglund et al, 2007). Many of these signals involve the 
oncogenic proteins KRAS, BRAF, PIK3CA and the tumour 
suppressor PTEN which are all downstream effectors of the 
epidermal growth factor receptor (EGFR) (Siena et al, 2009). 
Treatment targeting EGFR has been found to be efficient only if no 
mutations are found in KRAS or BRAF (Lievre et al, 2006). Still 
all patients with wild-type KRAS and BRAF do not respond to 
treatment (Amado et al, 2008; Bardelli and Siena, 2010; Tol et al, 
2010). PIK3CA and PTEN have been suggested to harbour 
aberrations in 30-40% of all sporadic CRC cases (Samuels and 
Ericson, 2006; Frattini et al, 2007), which might explain part of this 
resistance. A recent study suggested that mutations in PI3K 
catalytic subunit (PIK3CA) may carry prognostic information in 
tumour stage I-III (Ogino et al, 2009), and that PIOCA/PTEN 
deregulation, in addition to KRAS and BRAF mutations, may be a 
biomarker of resistance (Perrone et al, 2009; Sartore-Bianchi et al, 
2009). Consequently, Sartore-Bianchi et al (2009) introduced the 
Quadruple index as a factor taking aberrations in these four factors 
into simultaneous consideration. Even though many studies are 
focusing on the molecules downstream EGFR to estimate 
benefit from EGFR blocking therapy, it is still not known how 
the mutations affect patient prognosis and tumour aggressiveness 
per se. 

Therefore, we have in the present study analysed the mutational 
status of KRAS, BRAF, PIK3CA and PTEN expression separately, 
and combined as Quadruple index, and correlated the results to 
patient survival. Additionally, we related mutation status to 
established molecular tumour characteristics such as MSI screening 
status and CIMP status. 



MATERIAL AND METHODS 



Patient selection. Colorectal cancer cases from two separate 
patient groups were included in the present study. Archival 
paraffin- embedded CRC tissue samples from a total 414 patients 
were included from the Colorectal Cancer in Umea Study 
(CRUMS), all collected during primary tumour surgery over 
the period 1995-2003 at Umea University Hospital, Sweden. 
All routinely stained sections were reviewed by one observer, 
who performed all histopathological classifications including stage 
and tumour type (mucinous or non-mucinous). Tissue blocks from 
the primary tumour were chosen for DNA extraction. When 
necessary the proportion of tumour cells was maximised by 
macrodissection and necrotic areas were avoided. Clinical data 
were obtained by reviewing the patient records and survival 
data were collected from the Swedish population registry during 
autumn 2012 with a median follow-up time of 113 months for 
patients still alive at the end of follow-up. 

From the Northern Sweden Health Disease Study (NSHDS), 
archival paraffin- embedded CRC tissue from a total of 197 patients 
was included. The NSHDS cohort consists of three separate 
cohorts: the Vasterbotten Intervention Project (VIP), the Northern 
Sweden WHO Monitoring of Trends and Cardiovascular Disease 
Study (MONICA) and the local Mammography Screening Project 
(MSP) (Hallmans et al, 2003). The CRC cases in the NSHDS 
cohort, protocols and selection principles used in the present study 
have previously been described in detail (Van Guelpen et al, 2006). 
Brief summary of subjects included in the NSHDS cohort: consists 
of both men and women in the age of 40, 50 and 60 years in VIP; 
both men and women ages 25-74 years in MONICA; and only 
women ages ~ 50-70 years in MSP. Within these cohorts, a total of 
226 CRC cases were identified and selected for a previous 
nested case- referent study (Van Guelpen et al, 2006). After 
exclusion of insufficient or unavailable tumour tissue samples, 



197 patients were available for mutation analysis in the NSHDS 
cohort. 

NSHDS patients were followed up until January 2008 with a 
median follow-up time of 102 months for patients still alive at 
the end of follow-up. Cancer- specific survival was collected from 
the Swedish population registry and patient records. Patients 
originally included in both cohorts were excluded from the 
CRUMS cohort and only reported once. 

The handling of tissue samples and patient data in this study has 
been approved by the local ethics committee of Umea University, 
Umea, Sweden. 

Mutational analysis of KRAS and PIK3CA exon 20. PCR 

conditions for KRAS: 50 ng DNA, 0.5 fig primer, 10 mM 
dNTP, 1 mM MgCl 2 and 0.4U JumpStart Taq (Sigma, Stockholm, 
Sweden) in a total volume of 20 fil. PCR were run at 95 °C 10 min, 
95 °C 15 s, 65-55 °C ( - 1 °C/cycle) 72 °C 30 s (touchdown for 10 
cycles); 95 °C 15 s, 55 °C 15 s, 72 °C 30 s for 35 cycles and 72 °C 
10 min. Primers used: 

forward: 5' -tgtaaaacgacggccagtgagtttgtattaaaaggtactgg-3'. 

reverse: S'-caggaaacagctatgacctctgtatcaaagaatggtcct-S'. 

PCR conditions for PIK3CA exon 20: 50 ng DNA, 0.5 /ig primer, 
10 mM dNTP, 3mM MgCl 2 and 0.4U JumpStart Taq (Sigma, 
Stockholm, Sweden) in a total volume of 20 /A. PCR were run at 
95 °C 10 min, 95 °C 21s, 59 °C 21s, 72 °C 30 s for 40 cycles and 
72 °C 10 min. Primers used: 

forward: S'-tgtaaaacgacggccagtctcaatgatgcttggctctg-S'. 

reverse: S'-caggaaacagctatgaccatgctgttcatggattgtgc-S'. 

All primers were M13-tagged (forward: S'-tgtaaaacgacggccagt-S'; 
reverse: S'-caggaaacagctatg-S') to receive a more specific 
PCR product during the sequencing reaction. Sequencing 
was performed using Big Dye v. 3.1 according to the manufacture 
protocol, analysed in a 3730 xl DNA Analyser (Applied 
Biosystems, Stockholm, Sweden). The results were evaluated in 
SeqScape v2 1.1 (Applied Biosystem). 

BRAF V600E mutational analysis. Detection of BRAF V600E 
mutation was done with the Taqman allelic discrimination assay 
(reagents from Applied Biosystems), which has been described in 
detail elsewhere (Benlloch et al, 2006). 

Immunohistochemical analysis of PTEN expression. Specimens 
were fixed in 4% formaldehyde and embedded in paraffin, 
according to routine procedures at the Department of Clinical 
Pathology, Umea University Hospital, Sweden. Four micrometre 
sections were deparaffinized and rehydrated. Antigen retrieval 
treatment was executed using Borg solution (Biocare Medical, 
Concord, CA, USA) in a pressure cooker (2100 retriever, Biocare 
Medical). Primary monoclonal mouse PTEN antibody (Dako, 
Stockholm, Sweden, clone 6H 2.1, diluted 1 : 50) was used in a 
semiautomatic staining machine (intelliPATH FLX, Biocare Medical). 

The samples were evaluated for cytoplasmic staining, and were 
graded Oas no staining, las weak staining, and 2as moderate-strong 
staining. Loss of PTEN expression (graded as 0) was considered as 
abnormal while grade 1 and 2 was considered normal. Nerve tissue 
and blood vessels were used as positive internal controls in each 
sample. Cases without internal positive control staining were 
considered uninformative. 

A Quadruple index was created according to Sartore-Bianchi 
et al (2009), where negative specify cases where all selected genes 
(KRAS, BRAF and PIK3CA) were wild-type and normal expression 
of PTEN was seen. Quadruple index positivity indicates cases 
where at least one of the KRAS, BRAF or PIK3CA genes was 
mutated and/or loss of PTEN expression was found. 

Microsatellite instability screening status and CIMP status. 

Immunohistochemical analyses of mismatch repair proteins were 
performed as previously described (Dahlin et al, 2010). Briefly, 
expression of four mismatch repair proteins, MLH1, MSH2, MSH6 
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Figure 1 . The interrelationship between cases with mutations in KRAS, 
BRAF, PIK3CA and loss of PTEN expression in the NSHDS and the 
CRUMS cohorts. Total number of aberrations in NSHDS (A); KRAS 
(N = 30), BRAF (N = 31), PIK3CA (N = 3), PTEN (N=18); CRUMS (B); 
KRAS (N = 77), BRAF (N = 50), PIK3CA (N = 8), PTEN (N = 57). 
Patients with missing value in any of the marker were excluded 
from the Figure. 

and PMS2 were analysed in formalin-fixed and paraffin -embedded 
human CRC tissue. Tissue samples lacking nuclear staining in 
tumour cells for at least one of these proteins were considered to 
have a positive MSI screening status, referred to as MSI. Negative 
MSI screening status based on immunohistochemical staining 
is referred to as MSS. 

Methylation analysis to determine tumour CIMP status 
was performed by the MethyLight method, with primer and 
probe sequences as previously described (Weisenberger et al, 2006; 
Dahlin et al, 2010). The per cent of methylated refence (PMR) 
value was caluculated for the eight genes included in the CIMP 
panel (CDKN2A, MLH1, CACNA1G, NEUROG1, RUNX3, SOCS1, 
IGF2 and CRABP1) (Dahlin et al, 2010), and a gene was considered 
positive for methylation when the PMR > 10 (Weisenberger 
et al, 2006). 

Tumours with no promoter hypermethylation were classified as 
CIMP-negative, 1-5 genes methylated as CIMP-low, and 6-8 genes 
as CIMP-high (Dahlin et al, 2010). 

Statistical analysis. Clinico -pathological characteristics were com- 
pared using Kruskal-Wallis tests for continuous variables and j 1 - 
tests, or Fisher's exact tests when observed or expected frequencies 
were less than five for categorical variables. For cancer- specific 
survival analyses, Kaplan-Meier plots were used, and differences 
between groups were tested by log-rank tests. Cancer-specific 
events were defined as death with known disseminated or recurrent 
disease, and cases were censored at the end of follow-up or at time 
of death by other causes. 

Patients in CRUMS who were deceased with postoperative 
complications within 1 month after surgery (n = 16) were excluded 
from the survival analyses. Deaths due to postoperative 



complications were not recorded in NSHDS, but only four patients 
died within 1 month of surgery. To take into consideration other 
clinico-pathological factors, multivariate Cox proportional hazard 
models were used. For multivariate analyses, we analysed 
Quadruple index, KRAS and BRAF and not PIK3CA and PTEN, 
as the latter two were not significantly associated with prognosis in 
univariate analyses. The adjusting variables were selected if they 
affected the risk estimates for KRAS and BRAF > 10% in bivariate 
analyses. The final multivariate model included sex, age at 
diagnosis, stage and tumour site. Other factors tested, but not 
meeting the criteria for inclusion in the multivariate analyses were 
aberrant p53 protein expression, mucinous histologic tumour type, 
preoperative radiotherapy and adjuvant chemotherapy. Micro- 
satellite instability screening status and CIMP status were also 
tested but excluded due to small subgroups and thereby loss of 
statistical power. All statistical tests were conducted using PASW 
Statistics 18 (SPSS Inc., Chicago, IL, USA). 



RESULTS 



Quadruple index in relation to clinico-pathological variables. 

We analysed each mutation (KRAS, BRAF and PIK3CA) and PTEN 
expression as well as the Quadruple Index, in tumours from 197 
patients in the NSHDS and 414 patients in the CRUMS cohort 
with respect to different clinico-pathological characteristics 
(Tables 1A and IB). Seven different activating mutations in codon 
12 and 13 were analysed in KRAS, and the mutation frequency was 
17.9% in the NSHDS and 19.5% in the CRUMS cohort. BRAF was 
observed in 17.9 and 13.2% in each study population respectively 
(Tables 1A and IB). When combining results from the four studied 
factors, only two patients had both BRAF and KRAS mutated in the 
NSHDS cohort (Figure 1A), while BRAF and KRAS mutations were 
mutually exclusive (Figure IB) in the CRUMS cohort. Four 
different mutations were analysed in PIK3CA, exon 20, where the 
mutation frequency was 2.2% in both cohorts. Loss of PTEN 
expression was found in 12.5% in the NSHDS and 14.1% in the 
CRUMS cohort (Tables 1A and IB). In the NSHDS cohort mutated 
KRAS and BRAF tumours were associated with right colon 
location, most distinct for BRAF (NSHDS; P< 0.001). In the 
CRUMS cohort, BRAF mutant tumours were significantly 
correlated to older age (CRUMS; P = 0.017) and right colon 
location (CRUMS; P< 0.001), while KRAS mutations were 
significantly associated with higher tumour stage (CRUMS; 
P = 0.030). BRAF mutations were most prevalent in mucinous 
tumours (Tables 1A and IB). 

The frequencies of Quadruple index positivity were 48.3% in the 
NSHDS and 44.0% in the CRUMS cohort. Quadruple index 
positivity was correlated significantly to right colon location in 
both patient groups (NSHDS and CRUMS; both P< 0.001). 
Quadruple index positivity, BRAF mutations and loss of PTEN 
expression were significantly associated with higher tumour stage 
in the NSHDS, but not in the CRUMS cohort (Tables 1A and IB). 

Quadruple index in relation to MSI screening status and CIMP 
status. Tables 2A and 2B shows Quadruple index and each 
mutation (KRAS, BRAF and PIK3CA) and PTEN expression in 
relation to both MSI screening status and CIMP status in the 
NSHDS and the CRUMS cohort. Quadruple index positivity 
correlated significantly to CIMP-high status (NSHDS; P = 0.002 
and CRUMS; P< 0.001) in both the NSHDS and the CRUMS 
cohort, and to MSI (CRUMS; P< 0.001) in the CRUMS cohort. 
KRAS mutations were more often seen in patients with MSS 
(NSHDS; P = 0.031 and CRUMS; P = 0.002) and CIMP-low 
tumours (NSHDS; P = 0.046 and CRUMS; P = 0.001). BRAF 
mutations were significantly associated with MSI (NSHDS; 
P< 0.001 and CRUMS; P< 0.001) and CIMP-high (NSHDS; 
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Abbreviations: CIMP = CpG island methylator phenotype; MSS = micro-satellite stable; NSHDS = Northern Sweden Health Disease Study; MSI = microsatellite instability; Wt = wild-type. The 
following numbers of missing cases were present in NSHDS: CIMP status, 3; Quadruple Index, 25; KRAS mutation status, 18; BRAF mutation status, 1; PIK3CA mutation status, 11; PTEN mutation 
status, 13. Cases lacking nuclear staining of tumour cells for at least one of MLH1, MSH2, MSH6 or PMS2 were considered to have a positive MSI screening status (MSI). CIMP according to an 
eight-gene panel including CDKN2A, hMLHl, CACNA1G, NEUROG1, RUNX3, SOCS1, IGF2 and CRABP1; CIMP-negative, 0 genes hypermethylated; CIMP-low, 1-5 genes hypermethylated; 
CIMP-high, 6-8 genes hypermethylated. Kruskall-Wallis test was used for continuous variables, # 2 -test or Fisher's exact test used for categorical variables. 



Table 2b. Molecular c 


haracteristics of colorectal cancer cases in the CRUMS cohort 








N 


MSI 


MSS 


P-va/ue 


CIMP-negative 


CIMP-low 


CIMP-high 


P-va/ue 


Frequency (%) 


414 


62 (15.5) 


338 (84.5) 




209 (50.6) 


155 (37.5) 


49 (11.9) 




Quadruple Index 








< 0.0001 








< 0.0001 


Negative 
Positive 


227 (56.0) 
178 (44.0) 


19 (31.7) 
41 (68.3) 


201 (60.5) 
131 (39.5) 




142 (69.3) 
63 (30.7) 


82 (54.3) 
69 (45.7) 


3 (6.3) 
45 (93.8) 




KRAS 
Wt 

Mutant 


331 (80.5) 
80 (19.5) 


59 (95.2) 
3 (4.8) 


263 (78.3) 
73 (21.7) 


0.002 


174 (83.7) 
34 (16.3) 


111 (72.1) 
43 (27.9) 


46 (93.9) 
3(6.1) 


0.001 


BRAF 
Wt 

Mutant 


356 (86.8) 
54 (13.2) 


27 (44.3) 
34 (55.7) 


317 (94.6) 
18 (5.4) 


<0.0001 


206 (99.0) 
2(1.0) 


143 (92.9) 
11 (7.1) 


7 (14.6) 
41 (85.4) 


< 0.0001 


PIK3CA Exon20 








0.013 








0.006 


Wt 

Mutant 


396 (97.8) 
9 (2.2) 


55 (93.2) 
4 (6.8) 


328 (98.5) 
5(1.5) 




204 (99.0) 
2 (1 .0) 


150 (98.0) 
3 (2.0) 


42 (91.3) 
4 (8.7) 




PTEN 








0.719 








0.729 


Normal 
Loss 


352 (85.9) 
58 (14.1) 


52 (83.9) 
10 (16.1) 


286 (85.6) 
48 (14.4) 




178 (85.6) 
30 (14.4) 


134 (87.6) 
19 (12.4) 


40 (83.3) 
8 (16.7) 




Abbreviations: CIMP = CpG island methylator phenotype; CRUMS = Colorectal Cancer in Umea Study; MSI = microsatellite instability; MSS = microsatellite stable; Wt = wild-type. The following 
numbers of missing cases were present in CRUMS: CIMP status, 1; Quadruple Index, 9; KRAS mutation status, 3; BRAF mutation status, 4; PIK3CA mutation status, 9; PTEN mutation status, 
4. Cases lacking nuclear staining of tumor cells for at least one of MLH1, MSH2, MSH6 or PMS2 were considered to have a positive MSI screening status (MSI). CIMP according to an eight-gene 
panel including CDKN2A, hMLHl, CACNA1G, NEUROG1, RUNX3, SOCS1, IGF2 and CRABP1; CIMP-negative, 0 genes hypermethylated; CIMP-low, 1-5 genes hypermethylated; CIMP-high, 
6-8 genes hypermethylated. Kruskall-Wallis test was used for continuous variables, / 2 -test or Fisher's exact test used for categorical variables. 



P< 0.001 and CRUMS; P< 0.001). Mutations in the PIK3CA gene 
significantly correlated to MSI (CRUMS; P = 0.013) and CIMP- 
high (CRUMS; P = 0.006) in the CRUMS cohort, but showed no 
statistical significance in the NSHDS cohort. Loss of PTEN 
expression did not show significant correlation to MSI screening 
status or CIMP status in any of the cohorts. 

Survival analysis. Cancer- specific survival analyses revealed that 
Quadruple index positive cases had a significantly worse prognosis 



compared with negative cases in the NSHDS cohort (Figure 2A; 
univariate HR 1.98, 95% CI: 1.25-3.13). However, the Quadruple 
index positive cases had only a slightly poorer, but not statistically 
significant, prognosis in the CRUMS cohort (Figure 2B; univariate 
HR 1.22, 95% CI: 0.88-1.69). 

When analysing each gene separately only BRAF mutations 
turned out to be of prognostic value in the NSHDS cohort 
(Figure 2E), a result that retained statistical significant also in a 
multivariate Cox proportional hazard model (Table 3A). 



2158 



www.bjcancer.com | DOI:10.1038/bjc.2013.212 



Role of KRAS, BRAF, PIK3CA and PTEN in CRC 



BRITISH JOURNAL OF CANCER 



NSHDS 



CRUMS 



Table 3a. Cox regression of colorectal cancer cases in the NSHDS cohort 
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Figure 2. Cancer-specific survival analyses with respect to the 
Quadruple index and the KRAS, BRAF, PIK3CA and loss of PTEN 
expression separately. 



In the CRUMS cohort, on the other hand, only KRAS mutations 
were of prognostic value (Figure 2D), and this was seen also in 
multivariate analyses (Table 3B). Neither PIK3CA mutations, nor 
loss of PTEN expression were of prognostic significance in any of 
the two cohorts when analysed separately (Figure 2G-J). 



N 


Univariate 


Multivariate 




HR (CI 95%) 


HR (CI 95%) 


Quadruple Index 


172 


1.978 (1.251-3.128) 


1.308 (0.787-2.174) 


KRAS 


179 


1.325 (0.773-2.271) 


0.798 (0.443-1 .438) 


BRAF 


196 


2.428 (1.490-3.956) 


1.998 (1.165-3.426) 


PIK3CA Exon20 


186 


0.657 (0.091-4.739) 


0.285 (0.038-2.141) 


PTEN 


184 


1.555 (0.859-2.816) 


1 .289 (0.699-2.376) 


Abbreviations: C 


= confidence interval; HR = hazarc 


ratio, NSHDS = Northern Sweden 


Health Disease Study. HR determined by Cox proportional hazard models, adjusted for sex, 


age, tumour site and tumour stage. 





Table 3b. Cox regression of colorectal cancer cases in the CRUMS cohort 



N 


Univariate 
HR (CI 95%) 


Multivariate 
HR (CI 95%) 


Quadruple Index 


372 


1.220 (0.881-1.689) 


1.157 (0.827-1.619) 


KRAS 


378 


1.761 (1.220-2.542) 


1.485 (1.023-2.155) 


BRAF 


377 


0.843 (0.508-1 .397) 


0.914 (0.529-1.576) 


PIK3CA Exon20 


372 


0.000 (0.000-1.408 E + 122) 


0.000 (0.000-1 .088E1 69) 


PTEN 


377 


0.870 (0.531-1 .426) 


0.862 (0.519-1.431) 


Abbreviations: CI = confidence interval; CRUMS = 
HR = hazard ratio HR determined by Cox proport 
age, tumour site and tumour stage. 


= Colorectal Cancer in Umea Study; 
onal hazard models, adjusted for sex, 



Survival analyses stratified for MSI screening status and CIMP 
status. Patients with Quadruple index positive tumours with 
MSS (NSHDS; P = 0.002), or CIMP-low (NSHDS; P = 0.022) or 
CIMP-high tumours (CRUMS; P = 0.042) had a worse prognosis 
than Quadruple index negative cases. Cancer- specific survival 
analyses stratified for KRAS and BRAF is shown in Figure 3. 
Patients with tumours harbouring BRAF mutations together 
with MSS (NSHDS; P= < 0.001) (Figure 3G) or CIMP-low 
(NSHDS; P< 0.001) (Figure 30) showed an impaired survival 
in the NSHDS cohort. In the CRUMS cohort, tumours 
with KRAS mutations accompanied with MSS (Figure 3F) 
(CRUMS; P = 0.042) or CIMP-negative (CRUMS; P = 0.010) or 
BRAF mutations in CIMP-high tumours (CRUMS; P = 0.001) 
(Figure 3T) showed a poorer patient prognosis. Owing to the loss 
of statistical power in these small subgroups, a multivariate model 
was not performed. 
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Figure 3. Cancer-specific survival analyses in the NSHDS and the CRUMS, stratified for KRAS or BRAF mutations, in relation to MSI screening 
status and CI MP status. 



DISCUSSION 



In this study archival CRC tissue from two different cohorts from 
Northern Sweden, NSHDS and CRUMS, were analysed regarding 



mutations in the genes KRAS, BRAF, PIK3CA and loss of PTEN 
expression. All four aberrations investigated in this study are part 
of the same signalling pathway, downstream the EGFR, and to get 
an increased understanding for how these factors are intercon- 
nected in CRC, a Quadruple index as suggested by Sartore-Bianchi 
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et al (2009) was created, where Quadruple index positive tumours 
had at least one mutation in any of the genes KRAS, BRAF, 
PIK3CA and/or loss of PTEN protein expression. 

We found a shorter cancer- specific survival in patients with 
Quadruple index positive tumours in the NSHDS cohort, but the 
Quadruple index was not statistically significant in the CRUMS 
cohort. Analysing each gene separately revealed that only 
mutations in the BRAF gene had a significant prognostic value 
in the NSHDS cohort, especially in combination with MSS or 
CIMP-low. Only KRAS mutations, on the other hand, indicated a 
significantly poorer patient prognosis in the CRUMS cohort, 
especially together with MSS or CIMP-negative tumours. Aberra- 
tions in PIK3CA and PTEN did not add significant prognostic 
information. Therefore, our results do not support the use of the 
full Quadruple index but instead emphasise the prognostic 
information in KRAS and BRAF mutation status. 

Taken together, these results indicate that the establishment of 
molecular subgroups of CRC based on KRAS and BRAF mutation 
status can supply important information, not only in prediction of 
the EGFR-treatment response but also in prediction of patient 
prognosis. Importantly, KRAS and BRAF mutations are nearly 
mutually exclusive in CRC (Jakubauskas and Griskevicius, 2010; Li 
et al, 2011; Krol et al, 2012). 

The finding of contrary significances for KRAS and BRAF 
mutations in the two cohorts is not easily explained. However, it 
should be noted that the composition and the underlying design of 
the two cohorts differs significantly. For example, NSHDS consists 
of more women than men as a direct result of including the 
Mammary Screening Project as one of the three subcohorts, 
and BRAF mutations have more often been reported in women 
(Ogino et al, 2012). Furthermore, the age distribution also differs 
between the two cohorts and might have impact on the results. 
Not only the KRAS and BRAF mutations, but also molecular 
characteristics such as MSI screening status and CIMP status, are 
well known to correlate with the age and sex distribution 
(Nosho et al, 2009; Kalady et al, 2012). The contradictory results, 
however, emphasise a need for further larger studies on this topic. 

One of the main strengths of this study was the two large, non- 
overlapping, patient groups, which were both from the same 
northern Swedish population but had different recruitment 
protocols, age range and sex distributions. The patients in the 
present study were generally diagnosed previous to the broad 
introduction of many novel therapies, including successful 
resection of liver metastases, into clinical practice. Treatment was 
thus fairly homogeneous within each tumour site and stage. 
Residual confounding effect due to differences in treatment is 
therefore unlikely. It is not possible, however, to analyse the 
predictive value of mutations with respect to EGFR-blocking 
therapy in our patient cohorts due to the lack of such treatment 
during the cohort recruitment. Instead, the two cohorts include all 
tumour stages and are suitable for studies on tumour aggressive- 
ness and prognosis. 

The present study is, to the best of our knowledge, the largest 
study today on this subject. Despite the use of two patient cohorts, 
a limitation is, however, still the relatively low number of patients, 
especially when analysing somewhat rare subgroups (e.g., PIK3CA 
mutations, MSI cases or CIMP-high cases). The fact that we could 
not detect any correlation between loss of PTEN expression or 
PIK3CA mutations and patient prognosis makes us speculate that 
the need for analysing all four genes, as in the Quadruple index, 
might be unnecessary when prognosticating cancer-specific 
survival. There are, however, contradictory reports indicating that 
both PIK3CA mutations and loss of PTEN protein expression do 
affect patient prognosis (Sawai et al, 2008; Li et al, 2009; Jang et al, 
2010; Liao et al, 2012). 

The mutation frequencies of each analysed gene found in this 
study were in general similar to previous reports (Rako et al, 2012; 



Soeda et al, 2012), except for the KRAS gene. We report a 
frequency of about 20%, while several other reports have reported 
frequencies of 30-40% (Kim et al, 2012). The low mutation 
frequency of KRAS in our studied populations can have several 
explanations. Our patient cohorts have a rather high proportion of 
rectal cancers, and rectal cancers have a lower KRAS mutation 
frequency than colon cancers. Technical differences between 
studies are another likely explanation, and here we have not 
analysed KRAS mutations in exon 61. Furthermore, most studies 
reporting the frequency of KRAS mutations have studied only 
metastatic CRCs, and KRAS-mutated CRC might be more 
aggressive than their wild-type counterparts. 

Previous reports on PIK3CA mutation frequencies in CRC have 
varied considerably. In this study we report a frequency of about 
2%. However, we have only analysed mutations in exon 20 in 
PIK3CA, not exon 9, based on recently published data showing that 
only mutations in exon 20 have a prognostic value (De Roock et al, 
2010; Farina Sarasqueta et al, 2011), probably as this 
exon translates the kinase domain of PIK3CA. Additionally 
Muller et al (2007), recently found a PIK3CA pseudogene spanning 
exons 9-13 located on chromosome 22, which might be the reason 
for such a high reported frequency of PIK3CA exon 9 mutations. 

In conclusion, by the use of two patient cohorts we show that 
mutations in the KRAS and BRAF genes are of prognostic 
importance in colorectal cancer. However, adding information 
on mutation status of PIK3CA and loss of PTEN does not add 
significant prognostic information. These results suggest that 
establishment of molecular subgroups based on KRAS and BRAF 
mutation status is important and should be considered in future 
prognostic studies in CRC. 
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